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We review the critical behavior of nonequilibrium systems, such as directed percolation 
(DP) and branching-annihilating random walks (BARW), which possess phase transi- 
tions into absorbing states. After reviewing the bulk scaling behavior of these models, 
we devote the main part of this review to analyzing the impact of walls on their critical 
behavior. We discuss the possible boundary universality classes for the DP and BARW 
models, which can be described by a general scaling theory which allows for two inde- 
pendent surface exponents in addition to the bulk critical exponents. Above the upper 
critical dimension dc, we review the use of mean field theories, whereas in the regime 
d < dc, where fluctuations are important, we examine the application of field theoretic 
methods. Of particular interest is the situation in d = 1, which has been extensively in- 
vestigated using numerical simulations and series expansions. Although DP and BARW 
fit into the same scaling theory, they can still show very different surface behavior: for 
DP some exponents are degenerate, a property not shared with the BARW model. More- 
over, a "hidden" duality symmetry of BARW in d = 1 is broken by the boundary and 
this relates exponents and boundary conditions in an intricate way. 

Keyuiords: Phase transitions; Nonequilibrium systems; Absorbing states; Directed per- 
colation; Branching-annihilating random walks; Surface and bulk critical behavior. 
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1. Introduction 

As is well known from the study of equilibrium statistical mechanics, boundaries 
have important effects on systems close to criticalityE'ErH Quantities measured close 
to a wall can scale differently than in the bulk where, for a given bulk universality 
class, various boundary universality classes are possible. An analysis of the effects 
of boundaries is clearly very important if one wishes to apply the theory of critical 
phenomena to real physical systems. Although the theory of equilibrium surface 
critical phenomena is now well understood, the equivalent theory for nonequilibrium 
systems is in a much less advanced state. It is the purpose of this review to describe 
recent progress in extending the theory of surface critical phenomena to certain 
nonequilibrium problems. 

The most prominent example of such a nonequilibrium dynamic system is di- 
rected percolation (DP), which is the generic model for systems with a continuous 
phase transition from an active to an absorbing state from which the system cannot 
escape. DP describes the directed growth of a cluster (i.e. growth in a preferred 
spatial direction or along the time axis), where the growth rate is governed by a 
microscopic growth probability p. For probabilities below a critical value, p < pc, 
the cluster will die after a finite time, which means that the system gets trapped in 
the vacuum — the unique empty state. On the other hand, for high enough growth 
probabilities p > Pc, there is a finite probability that the cluster will always remain 
active. Exactly at ii = Pc, there is a critical phase transition from the active into 
the absorbing statcu A whole range of systems possessing a phase transition from a 
non-trivial active phase into a unique absorbing state fall into this universality class. 
Some examples include epidemics, chemical reactions, interface pinningAiepinning, 
the contact process, polynuclear growth, and certain cellular automataJju 

Although most models with absorbing states fall into the ubiquitous DP univer- 
sality class, there are important exceptions. For instance, the model of branching- 
annihilating random walks with an even number of offspring (BARW) exhibits quite 
different behavior0ij§ and defines a separate universality class. Other models in this 
class (at least in one spatial dimpnsinn d = 1) include certain probabilistic celluLai 
automataJi3 monomcr-dimcr models ,PrMr'i non-cquiliiii'ium kinetic Ising modelsjlil 
and generalized DP with two absorbing states (DP2).lla These models escape from 
the DP universality class by possessing an extra conservation law or symmetry: for 
the BARW model, a "parity" conservation of the total number of particles modulo 
2; for the other models, an underlying symmetry between their absorbing states. 
Models with an infinite number of absorbing states are believed to belong to the 
BARW class if the absorbing states can be arranged into two symmetric groups 
(see Ref. and references therein). On the other hand, if no higher symmetries 
between the absorbing states are present, then such models will belong to the DP 
class. 

In the present review we focus our attention on the impact of walls on systems 
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with nonequilibrium critical phase transitions into absorbing states.EHlBllilli'EliijiiSiijii'BiBii 

We wiU describe the boundary critical behavior of DP and BARW by using a variety 
of methods, including mean field theory, scaling theory, field theory, exact calcu- 
lations, Monte-Carlo simulations, and series expansions. The rest of this review is 
organized as follows: In Section ^ we introduce percolation and reaction-diffusion 
models for DP and BARW for bulk systems (without boundaries). Then in Section || 
we quickly review the general bulk scaling behavior of such systems. In Section ^ we 
briefly discuss the surface critical behavior of equilibrium systems, with a particular 
emphasis on the semi-infinite Ising model. In Section ^ we then introduce walls into 
our nonequilibrium models and discuss in detail some possible boundary conditions 
in rf = 1. In Sections ^, ^ ||, and |9| we review the various theoretical techniques 
used to analyze the surface critical behavior of our nonequilibrium models, includ- 
ing mean field, scaling, and field theories and also exact calculations. In Section 10 
we present results for critical exponents based on numerical work. In Section 11 
we review recent work in other directions concerning surfaces in DP-like systems. 
Finally we give a summary and outlook in Section ^ 

2. Bulk DP and BARW 

In this section we introduce the models for DP and BARW. We describe both 
reaction-diffusion versions of these models as well as the probabilistic cellular au- 
tomata frequently used in simulations in = 1. 



2.1. DP 

The easiest way of introducing DP is to consider bond percolation on a directed 
lattice. The update rules for bond DP in d spatial dimensions are then easily 
defined: for each site at time i, form bonds with probability p to the neighboring 
sites at time f -I- iD 

As is welllmawn, various reaction-diffusion models also fall into the DP univer- 
sality class.t3E3H The simplest of these is defined by the following reaction scheme: 

A ^ A + A with rate a 

A + A^A with rate A (1) 
A ^ (d with rate n, 

where the identical particles A otherwise perform simple random walks with diffu- 
sion constant D. 

In d = 1, the DP imi versa.l ity class has also frequently been studied by using 
the Domany-Kinzel model.aS In this cellular automata model each site can either 
be active or inactive and the probability for site i to be updated to state Si,t+i at 
time t -1-1 is given by an update probability P{si,t+i\si-i^t, Si+i,t)- See Figure |l| for 
a typical lattice configuration and Figure ^ for the update rules. An example of a 
cluster grown from a single seed according to these rules is shown in Figure Sa. 
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Fig. 1. Directed Pereolation in terms of the Domany-Kinzel model, where time flows vertically 
downwards. Black sites are active (A) and white ones inactive (/). The state of each site at time 
t + 1 depends on the states of the neighboring sites at time t. 



„o 

1 



Fig. 2. Update probabilities for DP in terms of the parameters < p, g < 1, where we have q = 
p{2 — p) for bond DP and q = p for site DP, respectively. Probabilities for the other configurations 
follow from left-right symmetry and from P(A | . . .) + P{I | . . .) = 1. 




(b) 






Fig. 3. (a) A DP cluster and (b) a DP2 cluster both grown from a single seed in the bulk. 
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2.2. BARW 

The BARW mfldel is defined by a reaction-diffusion system, with the foliowing 
reaction processesQ-Ha 

A (m + l)A with rate a^^ 
A + A^d) with rate A, (2) 

where the identical particles A otherwise perform simple random walks with dif- 
fusion constant D. For m odd, the above model is known to belong to the DP 
universality class, however for m even we have what we refer to as the BARW uni- 
versality class. Unless otherwise specified when we refer to the BARW model we 
will be referring to the even m case. 

The BARW class has also been studied in d = 1 using a generalized Domany- 
Kinzel model (sometimes called DPttO'EII) introduced by Hinrichsen.ta In this model 
each site can be either active or in one of n inactive states. For n = 1 the update 
rules are identical to those of the Domany-Kinzel model in Figure ^ but for n > 2, 
the distinction between regions of different inactive states is preserved by demanding 
that they are separated by active ones. An example of a DP2 cluster is shown in 
Figure The DP2 model has two symmetric absorbing states in which the system 
can become trapped. The update probabilities for d = 1, where DP2 is known to 
belong to the BARW universality class,E3 are given in Figure ^. 
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Fig. 4. Update probabilities for DP2: black sites are active {A), whereas white and grey sites are 
in the inactive states Ii and I2, respectively. Probabilities for the other configurations follow from 
left-right symmetry and from P{A | . . .) -|- P(/i | . . .) -|- P{l2 | . . .) = 1. 
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3. Bulk Scaling 

The growth of both DP and BARW clusters in the bulk close to criticality can be 
summarized by a set of independent exponents. A natural choice is to consider iy± 
and v\i which describe the divergence of the correlation lengths in space, 

a-\Ar\ (3) 

and time 

ClI-lAp''". (4) 

Here the parameter A describes the deviation from the critical point. For DP and 
DPn, A = pc — p, whereas for the reaction-diffusion models in mean field theory 
A = /i — cr for DP, and A = —mam for BARW. We also need the order parameter 
exponent (3, which can be defined in two a priori different ways: it is either governed 
by the percolation probability (the probability that a cluster grown from a finite 
seed never dies), 

P(< ^ oo,A) - |A|^-^ A<0, (5) 
or by the coarse-grained density of active sites in the steady state, 

n(A) - |A|''<^"% A<0. (6) 

When A < the system is said to be in an active state, whereas for A = 
the system is critical (with an algebraically decaying density), and for A > (if 
applicable) the system is either inactive (DP) or again critical (BARW).0 For the 
case of DP, it is known that [3 is unique: /3sccd = /3dcns in any dimension, both 
above and below ^e upper critical dimension dc =r^_This follows from time- 



reversal symmetrynd and field theoretic considerationsE§E3 and has been verified by 
extensive numerical work. The relation also holds for BARW in c? = 1, a result 
first suggested by numerics and now backed up by an exact duality mappingHa 
However, this exponent equality is certainly not always true: if we consider the 
BARW mean field regime valid for spatial dimensions d > dc — 2, then the system 
is in a critical inactive state only for a zero branching rate, where the density 
decays away as a power law. However, any nonzero branching rate results in an 
active state, with a nonzero steady state densityla (see Figure^). As the branching 
rate is reduced towards zero, this density (^) approaches zero continuously with the 
mean field exponent /3dcns = 1- Nevertheless, for c? > 2, the survival probability 
(H) of a particle cluster will be finite for any value of the branching rate, implying 
that /3seed = in mean-field theory. This result follows from the non-recurrence of 
random walks in d > 2. 

Field theoreticalhi DP is believed to be satisfactorily understood — the appro- 
priate field theoryEj'Eil (sometimes called Reggeon Field Theory, see Section 8.1) is 



well under control and the exponents have been computed to two loop order in an 
e = 4 — d expansion.E3 However, for the case of BARW (see Section |8.2|), a descrip- 
tion of the d = 1 case poses considerable difficulties for the field theory.Ll These stem 
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from the presence of two critical dimensions: dc — 2 (above which mean-field theory 
applies) and d'^ w 4/3. For d > d'^ the behavior of Figure ^ holds, i.e. an active 
state results for any nonzero value of the branching (t,„ , whereas for d < d'^ the sys- 
tem is only active for (t,„ > CTm, critical ,Q as shown in Figure This means that the 
physical spatial dimension d = \ cannot be accessed using controlled perturbative 
expansions down from the upper critical dimension dc = 2. Furthermore, for the 
Um < cTjn, critical region, the system is not inactive (in the sense of an exponentially 
decaying density). Instead this entire phase is controlled by the annihilation fixed 
point of the A + A ^ % process, where the density decays away as a power law. 
Hence this phase should rather be considered as still being critical. 




Fig. 5. Schematic bulk behavior for BARW of the density n as a function of the branching rate 
(Tm for (a) d > 2 and (b) d = 1. 

Despite the problems associated with BARW for d < d'^,we can still put forward 
a general scaling theory for DP and BARW, valid both above and below their critical 
dimensions. However, we must retain a possible distinction between /3sccd and /Sdons- 
For example, the average lifetime (t) of finite clusters can be derived from the scaling 
form for the survival probability 

P(t,A) = |A|''=-V(VC||)- (7) 

We then find 

{t) = J tP{t,A)dt^\A\-\ (8) 

where 

T = iy\\ - /3socd- (9) 

The appropriate scaling form for the density n(x,i), given that the cluster was 
started at x = 0, i = 0, is 

n(x,t. A) = |A|Aocd+/3de„=j(^/^^^i/^||). (10) 

Notice that rotational symmetry about the seeding point x = implies that the 
spatial coordinates enter the scaling function only as a; = |x|, the distance from 
the seeding point. Using the expression (|l^) we see that the average mass of finite 
clusters scales as 

(s) = J ?2{x,t,A) d'^xdt'^ \ A\-"' (11) 
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where 7 is related to the other exponents via the foUowing hyperscahng relation: 



+ dl^± = /3socd + /3dcns + 7- 



(12) 



Note that Eq. ( p^ ) is consistent with the distinct upper critical dimensions for 
BARW and DP. Using the above mean field values for BARW and u± ~ 1/2, 
i^ll = 1, and 7 = 1, we verify c?c — 2. In contrast, for DP one has the mean-field 
exponents /3dcns — Accd = 1, 7 = 1, and dc — 4. In Tables and ^ we list the 
exponents for DP (see Refs. |l|,||,|| , ^ , ^0|j4l|j4^ , ^ , ^ and references therein) and 
BARW (see Refs. p^,p5|,p6| and references therein). 



Table 1. Critical exponents for DP. 





d = l 


d = 2 


d = 3 


Mean Field 


/^dcns Pseed 


0.276 486(8) 


0.583(4) 


0.805(10) 


1 




1.733 847(6) 


1.295(6) 


1.105(5) 


1 




1.096 854(4) 


0.733(4) 


0.581(5) 


1/2 


T 


1.457 362(14) 


0.711(7) 


0.298(12) 





7 


2.277 730(5) 


1.593(7) 


1.232(12) 


1 



Table 2. Critical exponents for BARW. 





d=l 


Mean Field 




0.922(5) 


1 




0.93(5) 







3.22(3) 


1 




1.84(2) 


1/2 


T 


2.30(3) 


1 


7 


3.22(3) 


1 



4. Equilibrium Surface Critical Behavior 

We will begin our discussion of boundaries in critical systems byJpriefly reviewing 
the theory of surface critical phenomena for equilibrium systems,u with particular 
emphasis on the Ising model. Defined on a half space, the semi-infinite Ising model 
serves as the canonical example of a system exhibiting surface critical behavior. 

Each bulk Ising spin has the same number 2d of nearest neighbors. However, 
by cutting the lattice in half, spins next to the wall possess fewer neighbors than 
those in the bulk and are therefore more weakly coupled. At high temperatures T, 
the system is disordered and thus correlations are short ranged. Hence the effects 
of the wall are localized to a thin layer along the wall and decay away exponentially 
into the bulk with a length scale governed by the bulk correlation length. However, 
as the temperature is lowered the system will become critical at T = Tc, the critical 
temperature of the infinite system. At this point, where the correlation length 
diverges, the decay of the perturbation introduced by the wall is algebraic, and 
hence the impact of the wall on the scaling behavior is very marked. The two- 
point spin-spin correlation function decays with a new exponent along the wall. 
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independent of the two exponents required to describe bulk Ising criticality. In this 
case, the critical behavior of the surface is referred to as the ordinary transition. 
In two dimensions, the boundary critical e xpo nent governing the decay of the two 
point correlation function is known exactlyJ£3 

We could also imagine replacing all the spin couplings at the wall by stronger 
ones than in the bulk. Starting from the disordered state and lowering T, we will 
now have a situation where the strongly coupled wall spins may potentially undergo 
a phase transition already at T > Tc- In this case, the wall spins order independently 
of the bulk which remains disordered at this temperature. Of course, if the boundary 
is of sufficiently low dimension (e.g. the one dimensional boundary of a semi-infinite 
2d Ising model) then such ordering cannot occur at finite temperatures. But, for 
= 3, the surface is two dimensional and in these circumstances can undergo a 
so called surface transition. This surface transition (if it exists) lies in the same 
universality class as a d — 1 dimensional bulk transition. By further reducing the 
temperature, the bulk will then order at T = Tc in the presence of an already 
ordered surface. In this case, the critical behavior of the boundary is referred to 
as the extraordinary transition. This is quite different to the ordinary transition 
discussed above, where the surface was compelled to order through its coupling to 
the ordered bulk. At the extraordinary transition, where T = Tc, the system is again 
critical, but the correlations close to the wall differ from the case of the ordinary 
transition, and are governed by another independent boundary critical exponent. 

In principle we may tune the wall couplings such that the boundary transition 
takes place at T — Tc. Another name for this point where the bulk and boundary 
transitions coincide is the special transition. It is a multicritical point and connects 
the lines of ordinary, extraordinary and surface critical points in the phase diagram 
of bulk and wall couplings. 

For d = 2 the surface and special transitions cannot take place, since the surface 
is one-dimensional. However, an alternative way of ordering the surface layer is to 
introduce a magnetic field that only couples to wall spins. The boundary phase 
transition at T = Tc in the presence of such a field is called a normal transition. It 
turns out that the exact mechanism responsible for ordering the wall is unimportant. 
and thus the normal transition is in fact equivalent to the extraordinary transitionJla 
Hence, in two dimensions, we only need to distinguish between two universality 
classes. These can be accessed using the following boundary conditions: a free 
(open) boundary will give an ordinary transition and a fixed boundary will give an 
extraordinary transition. For both cases, the corresponding exponents are exactly 
known. Even the four-point functions governing the universal cros&jDver from wall 
to bulk correlations are exactly known from conformal field theory.Ea 

Another interesting aspect of the Ising model is the Kramers- Wannier duality, 
which is a mapping of the Ising model from disorder to order and vice versa.E3 The 
self dual point defines Tc. However, by introducing a surface, the self-duality at 
T = Tc is broken and it is straightforward to show that under the duality open 
boundaries are mapped to fixed boundaries and vice versa. Hence, the duality now 
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maps the two boundary universality classes onto one another. This observation is 
useful to bear in mind when analyzing the BARW model, where we will see that 
a self-duality broken by the wall again relates boundary universality classes in a 
subtle way. 



5. DP and BARW with Walls 

We now turn our attention to an understanding of the boundary critical prop- 
erties of the nonequilibrium DP and BARW models. We will begin by detailing the 
boundary conditions of these models. 

For the reaction-diffusion version of DP, the modification is simple, we simply 
allow for the DP reactions to occur at (potentially) different rates on the surface 
as compared to their values in the bulk. By tuning these boundary reaction rates, 
we can access the various boundary universality classes, in a similar way as was 
achieved in the Ising model by tuning the surface spin couplings. 

For the case of BARW, the situation is a little more complicated. The basic idea 
is that on the surface we may include not only the usual branching and annihilation 
reactions (|^) but potentially also a parity symmetry breaking A — > reaction. 
Depending on whether or not the — > reaction is actually present, we may then 
expect different boundary universality classes according to whether the symmetry 
of the bulk is broken or respected at the surface. We note that a somewhat similar 
situation in an equilibrium system has recently been analyzed in Ref. pU 

For the case of the d = 1 cellular automata introduced in Section there are 
several obvious boundary conditions. The simplest is just to cut off the lattice. 
This is equivalent to introducing boundary sites which are forced to be in one of 
the inactive states. We will refer to this case as the inactive boundary condition 
(IBC), and, for DP2, we choose inactivity of type 1 to the left of the boundary, see 
Figure Apart from imposing the state of these sites within the wall, the sites at 
the wall and those in the bulk are updated by the rules in Figures ^ and ^ 

Next we consider the reflecting boundary condition (RBC) where the wall acts 
like a mirror so that the sites within the wall are always a mirror image of those next 
to the wall, see Figure 0. For DP2, one can see that there is a qualitative difference 
between the IBC and the RBC. For the latter, regions of type-2 inactivity can get 
trapped at the wall and the only way for these regions to disappear is to wait for the 
cluster to return. On the other hand, for the IBC, such regions are never trapped. 

A third type of boundary condition is the active boundary condition (ABC) 
where the sites within the wall are forced to be active, see Figure ||. In this case the 
cluster will never die completely as the wall will always be active and can always 
induce new clusters. 

Although it is not readily apparent, it turns out that these three boundary 
conditions exhaust all possible boundary universality classes for the BARW model 



in d = 1 . We will return later in Sections |6.3.2| and |6.3.4| to the question of which 
boundary universality class is to be associated with each of the IBC, RBC, and 
ABC boundary conditions. 
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Fig. 8. DP2 with an active boundary condition (ABC). 



6. Boundary Critical Behavior of DP and BARW 
6.1. Surface DP 

First of all, let us examine the effects of introducing a d — 1 dimensional wall at 
= [x = (x||,a;j_)] into a DP process. Note that the labels parallel (||) and 
perpendicular (_L) refer here to directions relative to the wall (and not relative to 
the time direction). An example of such a cluster grown close to a wall is shown in 
Figure |^. 




Fig. 9. DP clusters in the Domany-Kinzel model grown from a single seed (a) in the bulk and (b) 
next to an impenetrable wall. 
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Fig. 10. Schematic mean field phase diagram for boundary DP. The transitions are labeled by 
0=ordinary, E=extraordinary, S=surface, and Sp=special. 



Following on from our discussion of boundaries in equilibrium systems (see Sec- 
tion 0), it is not difficult to justify a schematic phase diagram for boundary DP (see 



Figure hd and also Rcf. hq). In Figure 10, A and As represent, respectively, the 



deviations of the bulk and surface from criticality. For A^ > and as A ^ 0, we see 
an ordinary transition, since in that case we expect the bulk to order in a situation 
where the boundary, if isolated, would be disordered. On the other hand, for A > 
and for Ag sufficiently negative, we expect the boundary to order even while the 
bulk is disordered, i.e. the surface transition. Then for A^ < and A ^ 0, the bulk 
will order in the presence of an already ordered boundary, i.e. an extraordinary 
transition for the boundary. Finally at A = A^ = 0, where all the critical lines 
meet, and where both the bulk and isolated surface would be critical, we expect a 
multicritical point, i.e. the special transition. 

The bulk exponents are, of course, unchanged by the presence of a surface and, 
furthermore, one can show that the correlation length exponents on the boundary 
are also the same as in the bulk.ElEJ In the following, for conciseness, we will concen- 
trate solely on the ordinary transition; more details on some of the other transitions 
can be found in Refs. 16 2^. At the ordinary transition, one finds just one extra 



independent exponent associated with the boundary: this can be taken to be the 
surface density exponent /3i_dGns- This is defined from the steady-state density at 
the wall, where we have 

n{xi_ = 0, A) - I Al'^i^don, ^ A < 0. (13) 

In principle one could also allow for a second type of surface f3i exponent, one 
defined from a survival probability for clusters started on the wall 

Pi(f ^ oo. A) - |A|''°-<=d, A < 0. (14) 

However, the surface exponents here show a similar pattern to their bulk coun- 
terparts and fulfill /^Psccd = /^Pdcns = ^ ^an be shown using a time-reversal 
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symmetry argument (of. Ref. |2^ ) or by a field-theoretic derivation of an appropriate 
correlation function (see Section and also Ref. |l9|). 

6.1.1. Mean field theory for surface DP 

In this section, we very briefly review some simple results for boundary DP at the 
mean field level, focusing again on the ordinary transition. A more general analysis 
can be found in Ref. (see also Appendix C of Ref. U) . The equation describing 
mean field DP with a surface is 

dtn = DV'^n -An- An^, (15) 

with the boundary condition 

Ddx^n\^^=o = A^nl^^^o- (16) 

Here the variable A = /i — cr is the difference between the rates for the A ^ and 
A ^ A + A processes. Similarly we have the surface variable Ag, and the bulk 
quadratic term is due to the reaction A + A ^ A. From the above Eq. ([T^), the 
bulk mean field exponents can easily be computed: = 1, i^j^ = 1/2, and (3 — 1. 
Furthermore, with the inclusion of a boundary, we see that the correlation length 
exponents are unchanged at the wall but the surface f3i exponents are altered. If 
we are interested in the mean field steady state, then we can replace Eq. (Ilq) with 



Dn" -An- W = 0, (17) 

where n" = cPn/dx\. The appropriate boundary condition ( p^ ) is given by 

D< = Asns, (18) 

where ng = n\x^=o, and n'^ — dn/dxj_\x^=o. Multiplying Eq. ( pT[ ) by n' and 
integrating, we have 

^Dn'^ - ^An^ - ^\n^ + C = 0, (19) 

where C is a constant of integration. Using the bulk results n' — 0, and n — (—A) / A 
for A < 0, or n = for A > 0, we have 

^--m - {¥s + [A < 0] (20) 

^ = -[i]"'ns{in, + ^f'" [A>0] (21) 

where we have also used the boundary condition (p^. Considering the ordinary 
transition where A^ > and A 0^, we expect n = |A|/A ^ ris, and thus 
Eq. ( po|) yields oc |Ap/^, giving the mean field value /J^ — 3/2. 



Directed Percolation and Systems with Absorbing States 1775 



6.2. Surface BARW 

For the case of BARW, we expect a similar picture to hold as that described for 



DP in Section 6.1. In particular, the expressions (|13|) and ( |14| ) for the steady state 
density and survival probabilities will also apply. However, unlike DP, we will see 
that the /3i,sGod and /?i,dons exponents are no longer equal. 

6.2.1. Mean field theory for surface BARW 

The surface phase diagram for the mean field theory of BARW (valid for d > dc = 2) 
is shown in Figure |ll|. Here a^n, o"™^ are the rates for the branching processes 
A — > (m + 1)^ in the bulk and at the surface, respectively, and /is is the rate for 
the surface spontaneous annihilation reaction yl — + 0. Otherwise, the labeling is the 
same as that for the DP phase diagram (see Figure ^0|). We summarize the main 
details of the phase diagram below; more details can be found in Ref. |2^. 

The first feature to note is that the bulk is either active ((t„j > 0) or critical 
(cm ~ 0), but never inactive. Hence, unlike DP, there is no possibility of finding a 
surface transition, where the surface is critical with the bulk inactive. For the case 
where am = /^s = 0, we expect that for any finite value of the surface branching, the 
surface will become active. This corresponds to the extraordinary transition with an 
active surface and critical bulk. On the other hand for (t,„ = and /i^. > 0, 

the density at an (isolated) surface would decay away exponentially quickly due to 
the A 9 reaction, and hence we have the ordinary transition. Consequently with 
cr„j — 0, but both /is and am^ non-zero, there should be a line of special transitions 
dividing the extraordinary and ordinary regions. This explains the general features 
of the phase diagram in Figure |ri|. 

At a more quantitative level, the mean field equation for BARW is very similar 
to that for DP: 

dtn = DV^n -An- An^ (22) 
with the boundary condition 

Ddj:^n\x^=o = Asn\x^=o- (23) 

However the values of the A, As parameters are now different: A = —mam and 
As = —mams~^tJ-s- For simplicity we will again focus only on the ordinary transition 
(details of the other transitions can be found in Ref. |2^). Clearly we expect the 
same mean field exponent /3i,dcns as in DPEj However, the mean field behavior of the 
/3i,socd exponent is very different from the corresponding behavior in DP. Consider 
placing two particles next to the surface at i = 0. From the recurrence properties of 
random walks we see that, regardless of the reaction rates on the surface or in the 
bulk, there is a finite chance that the two particles will never meet again. Hence 
the survival probability is nonzero and thus /3psood = in mean field theory. 
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Fig. 11. Schematic mean field boundary phase diagram for BARW. See text for an explanation of 
the labeling. 



6.3. Beyond mean field theory 



6.3.1. Surface DP 



We expect that the phase diagram shown in Figure 10 is generally valid for surface 
DP close to the upper critical dimension dc = 4. However, as first pointed out in Ref. 
p7| , in d = 1, where the surface is just a zero dimensional point, the phase diagram 
may look rather different. In that case, for an inactive bulk, net particle production 
is only possible at one point. Furthermore, since particles will be constantly lost into 
the bulk, where they will decay away exponentially quickly, it may not be possible 
to form an active surface state. Of course this conclusion assumes the absence of 
a surface reaction 9 ~* A, which is the analog of a surface magnetic field in the 
Ising model. If this reaction is included then a normal transition obviously becomes 
possible. 

6.3.2. Boundary condition classification for d ^ I DP cellular automata 

Following the analysis in the previous section, we can now attempt to assign uni- 
versality classes for the IBC, RBC, ABC cellular automata boundary conditions 
introduced earlier in Section |[ The ABC condition obviously behaves as if there 
existed a surface reaction equivalent to — > A, and thus it belongs to the nor- 
mal transition universality class. Numerically, the IBC, RBC, belong to the same 
universality class ,E2I which we identify as the ordinary transition. 

6.3.3. Surface BARW 

Next we turn our attention to the d = 1 phase diagram for surface BARW shown in 
Figure |l^. The phase diagram looks quite different from its mean field analog due 
in part to the shift of the bulk critical point away from zero branching rate, but 
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(a) 



Sp=' 



(b) 



Sp 



O 



Fig. 12. Schematic surface phase diagrams for BARW in d = 1 for (a) am < o"m, critical i and (b) 
o"m = Tm, critical- See text for an explanation of the labehng. 



also due to the absence of any extraordinary transition. Again, in the absence of 
the surface reactions of the form — > ^, this is due to the fact that excess particle 
production (with a finite reaction rate) at a zero dimensional surface is simply not 
efficient enough to generate an active state, due to leakage into the critical bulk. 
We now divide the phase diagram into two sections, as in Figure [l^ . 

^rn ^ '-^m, critical ■ 

In this case, the bulk will be controlled by the annihilation process ^ + A ^ 0, 
and the branching process will be everywhere irrelevant. Consider first the case 
with /is = 0, where the surface reaction A — > is not permitted. In that case, 
since both the wall and bulk will be controlled by the critical A + A ^ reaction, 
we expect an analog of the special transition. However, this will not be a special 
transition in the usual sense, since in this region it will not be possible to obtain 
an active state either on the surface or in the bulk by small changes in the bulk 
and/or surface branching rates. To emphasize this point, we label this regime as 
Sp* in Figure This simpler and analytically tractable case has already been 
extensively analyzed in Ref. (see also Section |8.3| ). Next we consider the case 
fis > 0, where the reaction A ^ is permitted on the wall. Since the branching is 
irrelevant in this region of parameter space, the reaction A ^ would give rise to 
an exponentially decaying density at an isolated boundary, and hence an inactive 
state. Consequently we expect an analog of the ordinary transition. However, for 
the same reasons as described above, we will label this regime as O* in Figure [T^a. 

'^m — ^m, critical ■ 

Since we are now at the bulk critical point, which borders the bulk active phase, 
we expect rather different behavior to that described above. In this case it will 
be possible to obtain an active bulk/surface state by small changes in the reaction 
parameters. For the case where /i^ = 0, with both the bulk and surface critical, we 
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expect a standard special transition (labeled as Sp in Figure [T2p) . On the other 
hand, if fi^ > 0, then the parity symmetry of BARW is broken at the boundary 
and we expect a critical bulk and (if isolated) an inactive surface, i.e. a standard 
ordinary transition (labeled as O in Figure |l2|b) . 

6.3.4. Boundary condition classification for d — 1 BARW cellular automata 

We can now discuss the relations between the boundary conditions for the d — I 
DP2 cellular automata introduced in Section ^ and our above classification of the 
boundary universality classes for BARW in d — 1. The key feature is whether the 
symmetry in DP2 between the two absorbing states in the bulk is preserved at 
the surface. This symmetry is the analog of the parity symmetry conservation in 
the reaction-diffusion BARW model, which allows the DP universality class to be 
escaped. As we discussed above, breaking the parity symmetry on the boundary 
leads to the ordinary transition. The analog of that situation is provided by the RBC 
model which breaks the symmetry between the absorbing states on the boundary 
and hence also belongs to the ordinary transition. On the other hand, the IBC 
model respects the symmetry at the boundary, and thus belongs to the special 
transition. Furthermore, the ABC model once again clearly belongs to the normal 
transition universality class. Hence, as claimed earlier in Section |^, we see that 
by using the IBC, RBC, ABC classification all the previously discussed boundary 
BARW transitions in c? = 1 can be accessed. 



7. Scaling Theory 

In this section we review the scaling theory for the survival probabilities and corre- 
lation functions for boundary DP and BARW. For simplicity we again restrict our- 
selves to the ordinary transition; additional details can be found in Refs. p^j2l| , p7| . 

7.1. Surface DP 



As we mentioned in Section 5.1, the most important point to emphasize is that there 
is only one boundary /3p exponent for DP: /3P — fi^g^^^ — /3pdcns- The survival 
probability (the probability that the cluster is still alive at time t) has the form 

Pi(i,A)=A'^?Vi(t/e||), (24) 

where the scaling function ipi is constant for t » .0 The mean lifetime of finite 
clusters 

(t) ^ |A|-? (25) 
then follows from the lifetime distribution — dPi/diJl^ yielding the exponent 

=1/11-/3°. (26) 

However, for z^y < /3P, the leading contribution to (t) will be a constant, such that 
the above scaling relation breaks down and is replaced by ri =0. 
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For the coarse-grained bulk density ni at the point (x, t) of a cluster grown from 
a single seed located next to the wall at x = 0, t = 0, one can make the scaling 
ansatz 

ni(x,t, A) = A^?+^h {x/U, t/Cll) • (27) 
This ansatz may be properly justified using the field theory analysis reviewed in 



Section 8.1, however a more intuitive justification of the prefactors may be given as 

qO 



follows. The first factor of A'^i comes from the probability that an infinite cluster 
can be grown from the seed. The second factor of A'' arises from the bulk scaling 
of activity in the active state, i.e., the (conditional) probability that the point (x, 
t) belongs to the infinite cluster grown from the seed (see also Ref. In contrast, 
if the density is measured at the wall, then the appropriate ansatz reads 

nn{x,t,A) = A^^Vii (x/^, VCn) , (28) 

as we pick up a factor A^? rather than A'^ for the probability that (x, t) at the 
wall belongs to the infinite cluster grown from the seed. 

By integrating the cluster density (|2^) over space and time, we arrive at the 
average size of finite clusters grown from seeds on the wall, 

is) ^ \Ar?, (29) 

such that 

+djy_L + 7?- (30) 

Hence, the surface exponent 71 is related to the previously defined exponents via a 
scaling law that naturally generalizes the usual d+l dimensional bulk hyperscaling 
relation 

i^H+du± = 2/3 + -f. (31) 



The results of numerical simulations with a wall (see Section 10) are in very good 
agreement with the hyperscaling relation (|3^) . 

Besides integrating the density (|2^), we can also integrate the density on the 
wall (^8|) over the d — I dimensional wall and time. This integration yields the 
average (finite) cluster size on the wall, 

(^1,1) - |A|-^"S (32) 

where 

+(d-l)i.i = 2/3f +7°!. (33) 

However, in higher dimensions {d ~ 2 being a marginal case) this relation is not 
fulfilled as it would predict a negative Ji'i- For this case, Ji'i = 0, reflecting a 
constant contribution to Eq. (^^, cf. the comment after Eq. (^). 
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7.2. Edge DP 

We next turn to briefly review the case of DP clusters started on an edge. It has 
been known for some time that the presence of an edge introduces new exponents, 
independent of those associated with the bulk or with a surface (see Ref. ^ for a 
discussion in the context of equilibrium critical phenomena, or Ref. in the context 
of percolation). 

Consider a system, where we allow the wall to have an edge with an angle a 
at ~ x± ^ 0. Hence, the edge is simply the d — 2 dimensional intersection of 
two id — 1 dimensional walls. By placing the seed next to this edge, the boundary 
exponent /Jf is replaced by the edge exponent l32{a) (where of course /3P = /32'(7r)). 
Following the same arguments as before, we have the new scaling ansatz for the 
cluster density 

n2(r,i. A) = A^?+^h {r/^, , (34) 

where r is the radial coordinate in a system of spherical polar coordinates centered 
on ccii^^ — — Q. This ansatz applies for directions away from the edge and 
the walls. By replacing /3 with /3p or /^p, we get the corresponding results for the 
density along the wall or the edge, respectively. Moreover, in analogy with Eqs. (^ ) 
and ( ^ ) for seeds on a wall, we obtain the average (finite) size (s) ~ |A|~'''2 of 
clusters grown from a seed next to an edge, by integrating Eq. (|4|) over space and 
time. This yields the hyperscaling relation 

v\\+dy^=P^ +13 + ^^. (35) 

Quoting the results from Ref. |l^ for a mean field calculation of the edge exponents 

7^ = 1 - 7r/2a, = 1 + 7r/2a, (36) 

we see that Eq. (^5|) is satisfied at the upper critical dimension dc = 4. Numerical 
estimates for the exponent j^^ from Ref. ^ together with the mean field values, 
are listed in Table ^ (see Section |l^ for further details on simulation methods). 

Table 3. Numerical estimates for the /Jj* exponents for edge DP together with the mean field 
values from Eq. (|36|). Recall that /3^(7r) = /3p. The bulk and d = 1 wall estimates are listed for 
reference. 



Angle (a) 


7r/2 


37r/4 


71" 


57r/4 


bulk 


PY (d = l) 






0.73371(2) 




0.276486(8) 


(d = 2) 


1.6(1) 


1.23(7) 


1.07(5) 


0.98(5) 


0.583(4) 


(MF) 


2 


5/3 


3/2 


7/5 


1 



7.3. Surface BARW 

In this section we briefly review the analogous scaling theory for surface BARW. 
Again, for brevity, we will only consider the ordinary transition. When writing down 
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this scaling theory we must now bear in mind the important distinction between 
the /3i,dcns and /?i,sccd exponents. For a seed placed on the wall at x = 0, t = 0, the 
scaling form for the survival probability has the form 

Fi(t,A) = |A|^°=-$i(i/^ll). (37) 

It is then straightforward to compute the average lifetime of finite clusters, (t) ~ 
Aj^'^i, where = — Piseed^ J^^t ^ ^f DP. 

For the coarse-grained bulk particle density ni at the point (x, t) of a cluster 
grown from a single seed located next to the wall at x = 0, t = 0, one can make the 

ansatz 

ni(x,t, A) = |A|'3°-d+^^-5i {x/a, t/C\\) ■ (38) 



As was the case for DP, the A-prefactor in Eq. ( |38D comes from Eq. (|3^) for the 
probability that an infinite cluster can be grown from the seed, and from Eq. (^) 
for the (conditional) probability that the point (x, t) belongs to this cluster. If, 
instead, the density is measured at the wall, we have 

nn{x, t, A) = |A|ft°=-+ft%„.^^^ (^/^_^, t/^^^) , (39) 

as we pick up a factor |A|''i. rather than |A|''''""= from the probability that (x, 
t) at the wall belongs to the cluster. The average size of finite clusters 

{s,)^\Ar?, (40) 



follows from integrating the cluster density ( pSD over space and time, where the 
exponent 7P is related to the previously defined exponents via 

+ dl^l_ = + /Jdcns + 7P • (41 ) 

Analogously, by integrating the cluster wall density ( ]39| ) over the (d— l)-dimensional 
wall and time, we obtain the average size of finite clusters on the wall 

(^i,i)-|A|-^?S (42) 

where 

z.|| +{d- = (3^,,^ + /SO^cns + 7m- (43) 

Note that if the 7 exponents obtained from Eqs. (O) and (^) are negative, then 
they should be replaced by zero in Eqs. ( pO| ) and (p2^. 

8. Field Theory 

In the following sections we briefly review the available field theoretic results for 
bulk and boundaryElH DP and BARW. These techniques allow for the systematic 
inclusion of fluctuation effects, important below the upper critical dimension d(.. 
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8.1. DP field theory 

The field theory describing bulk DP is very well known. The action is given by 

where (/)(x, t) is the "density" field, and where 0(x, t) is the response field. Simple 
power counting reveals that the upper critical dimension is dc = 4 below which 
fluctuation effects become important. Renormalization of the theory is standard, 
involving field, mass (A), diffusion constant (D) and coupling constant (u) renor- 
malizations. The critical exponents can then be computed perturbatively using an 



e = 4 — c? expansion, givinj 



J30 



Note that, in this case, the exponents /?socd and /3dons can be shown to be equal. 
Technically, this follows from the fact that the density and response fields renor- 
malize identically (see also Ref. p4| ). 

In the action ( |44[ ) one can integrate out the response field (^(x, t) and arrive at 
a Langevin equation for the local density 0(x, t): 

= DV^cbi^, t) - A0(x, t) - iu0(x, tf + ry(x, t), (46) 

with 

(r;(x, t)) = 0, (rKx, t)77(x', t')) = u0(x, t)5\^ - ^')5{t - t'), (47) 

where ?7(x, t) is a Gaussian noise term. The multiplicative factor (/'(x, t) in the noise 
correlator reflects the fact that (f> = is the absorbing state. 

The use of field theories tQ_study boundary nonequilibrium phase transitions 
was initiated by Jansscn et al.ll3 They showed that the appropriate action for DP 
with a wall at a;^ = is given by 5 = S'buik + 'S'surfaco, where 

= A.4*., (48) 

with the definitions c/Js = 0(x|| , xj^ — 0,t) and (f>s = (/'(x||,x_l = 0,t). The sur- 
face term S'surfaco corresponds to the most relevant interaction consistent with the 
symmetries of the problem and which also respects the absorbing state condition. 
Simple power counting indicates that [As] ~ k, where k denotes an inverse length 
scale. The presence of the wall implies the boundary condition at x± = of 

Dd:,^(l)\^ = As(j)s- (49) 

Using this boundary condition, we see that a boundary term of the form (pgdx^ips, 
although marginal from power counting arguments, is actually redundant. 

Since Ag ^ k, its renormalized value can only flow to one of three possible fixed 
points: or ±oo. These possibilities correspond to the various types of boundary 
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phase transitions. The flow A, — oo (oo) corresponds to the extraordinary (or- 
dinary) transition, whereas the unstable fixed point at As — corresponds to the 
multicritical special transition. In the following we will concentrate solely on the 
ordinary transition, since that is the only case to have been studied numerically 
(see Ref. ^ for a field theoretic analysis of the special transition). 

The only extra renorrnalization required by the presence of the wall is a sur- 
face field renormalization.ll3 The presence of this extra renorrnalization at the sur- 
face leads naturally to the existence of just one independent surface exponent 
/3P. Once again, the exponents /^Pggcd ^'^^ /^Pdcns ^^"^ ^e shown to be equaill 
/^?seod — /^Pdcns — 01 ^ similar to the result found in the bulk. This follows from 
the fact that the surface density and response fields renormalize identically. Fur- 
thermore it can be shown that the correlation length exponents are everywhere 
unchanged by the wall (see Refs. The exponent can again be computed 

using an e = 4 — d expansion yieldingll3 

P? = /3?secd = /?i°de„s = ^ - 5 + Oie'). (50) 



From the field theory of Ref.Jlq, it is not hard to verify that Eq. (30) is the appro- 
priate generalization of Eq. (^l|) , relating /Jf to 

7? = ^ + | + 0(e^). (51) 



8.2. BARW field theory 

We now briefly review the bulk field theory for the BARW reaction-diffusion system. 
In order to properly include fluctuation effects for BARW, one must be careful to 
include processes generated by a combination of branching and annihilation. In 
other words in addition to the process A ^ {m + 1)A, the reactions A ^ [m —A) A, 
. . . , A iA need to be included. These considerations lead to the full actional 



m/2 



dt 



(v3(x,t)[at-i?vXx, t) 



1=1 



0-2; [1 



^(x,t)2']^(x,i)V'(x,i) (52) 
A[l - ^(x, i)2]^(x, tf) - V'(x, r) - no^(x, 0)1 , 



written in terms of the response field ip{x,t) and the "density" field ?/'(x, t). Here 
D is the diffusion constant, A the annihilation rate, and a2i the branching rate for 
the process A ^ {21 + 1)A. Note that the final two terms in Eq. (52) represent, 
respectively, a contribution due to the projection state (see Ref. ^5|), and the initial 
condition (an uncorrelated Poisson distribution with mean no). Notice also that 
(for even m) the action (p2[) is invariant under the "parity" transformation 



V'(x,i) ^ -^/'(x,t), '0(x,t) 



-^(x, t). 



(53) 
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This symmetry corresponds physically to particle conservation modulo 2 and it 
enables the system to escape the DP universality class. 

Simple power counting on the action in Eq. ( |5^ ) reveals that the upper critical 
dimension is dc = 2. Close to d^, the renormalization of the above action is quite 
straightforward (here we again quote the results from Ref. ^). At the annihilation 
fixed point the RG eigenvalue of the branching parameter can easily be computed. 
To one loop order one finds = 2 — TO(m+ l)e/2 + O(e^), where e = 2 — d. Hence 
we see that the lowest branching process is actually the most relevant. Therefore, 
close to 2 dimensions, where the branching remains relevant, we expect to find an 
active state for all nonzero values of the branching (as was the case for the BARW 
mean field theory reviewed earlier). 

However, inspection of the above most relevant RG eigenvalue j/o-^ shows that it 
eventually becomes negative at d = « 4/3. In that case we expect a major change 
in the behavior of the system, since the branching process will no longer be relevant 
at the annihilation fixed point. The critical transition point is then shifted with 
the active state only being present for values of the branching greater than some 
positive critical value (as indicated in Figure ||b). Consequently, we see that there 
is a second critical dimension d^ w 4/3 whose presence immediately rules out any 
possibility of accessing the non-trivial behavior expected in d = 1 via perturbative 
epsilon expansions down from d = 2. Instead cruder techniques (such as the loop 
expansion in fixed dimension) must be employed.Q 

We now review the effects of introducing a surface into the BARW universality 
classEH and further allow for the possibility of a symmetry-breaking A $ reaction 
to take place (with rate /i^), but only on the surface. 

8.3. fis^O BARW field theory 

In this case only the branching process is relevant on the surface. However, as in the 
bulk, one must still be careful to include the surface branching processes generated 
by a combination of branching and annihilation. This leads to a full surface action 
of the form 

^surface = J d'^-^X^^ ^ dt ^21^ (l - l^f ) | , (54) 

where i^s = ip{^\\,x± = 0,t) and ■0s = ij{^\\,x± = 0,t). Notice that the parity 
symmetry ( p3| ) is preserved for the /is = model at the wall, as well as in the bulk. 
This boundary action leads to boundary conditions very similar to those discussed 
earlier for DP, of the form Ddx^ip\x_L — ^sV's, where is the surface mass, equal 
to —mam^ in mean field theory. 

Power counting on the above action reveals that the surface branching rates 
(T2is all have naive dimension [cr2i^] ~ k^, where k denotes an inverse length scale. 
However, close to, but below d = 2, this scaling dimension will be renormalized 
downwards (this can be seen physically as a result of processes like A 3A A 
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rendering the branching process less efficient). As a result of this rcnormalization, 
we expect the lowest generated process (i.e. with ^ = 1 in Eq. (p^)) will become 
the most relevant (as it was in the bulk). Nevertheless, despite this downward 
rcnormalization, close enough to d = 2, the scaling dimension of the most relevant 
coupling will remain positive, and thus it will flow to oo for all nonzero starting 
values. This state of affairs corresponds to the extraordinary transition, where the 
surface is active while the bulk is critical. On the other hand, at bulk criticality and 
with (T2s = 0, we have a multicritical special transition point. Hence, for /i^ = 
and close to d = 2 the basic structure of the phase diagram is unchanged from 
mean field theory, although the values of the exponents will of course be altered by 
the fluctuations. In this case, after writing down and solving the appropriate RG 
equations (exactly along the lines of Refs. ^,^6h, one can derive scaling results for 
the density similar to those quoted in Section [7.3[ 

The situation in d = 1 is rather different, partly due to the shift of the bulk 
critical point away from am = 0. This ensures that the boundary and bulk transi- 
tions in d = 1 are inaccessible to controlled perturbative expansions. Nevertheless 
we still expect the scaling dimension of all the a2i^ to be negative in d = 1, fol- 
lowing the downwards trend in the renormalization mentioned above. In that case 
surface branching is then irrelevant in d = 1 leading to the Sp* and Sp special 
transitions discussed earlier. Hence, if the above scenario is correct, we do not ex- 
pect to see an extraordinary transition in d = 1 for any finite value of the surface 
branching, since the surface branching will always be irrelevant. This conclusion 
was investigated numerically in Ref. |2^, where no evidence of an active surface state 
for am < fm.criticai was found even for very high values of the surface branching 
parameter in a fermionic lattice model in d = 1. 

8.4. Hsj^Q BARW field theory 

In this case the reaction A — > is now possible, but only at sites on the wall. In- 
cluding surface processes generated by a combination of branching and annihilation 
(i.e. the processes A — > mA, A ^ {m ~ l)A, . . ., A 2A), the full surface action 
becomes 

/ dt (V ai^{l ~ i,i)i,si^s + ^is{i's - l)i^s] . (55) 



Note that the symmetry (^3[) is now broken by the surface term proportional to fig , 
which describes the A ^ reaction. Once again we have a boundary condition of 
the form Ddxj^tp\x± = AsV's, where now ~ —mams + f-s in mean field theory. 

The renormalization of the action ( |55| ) is now somewhat different from the jig ~ Q 
case. We again expect that we need only keep the lowest generated branching term 
on the surface, namely that with I — I'm Eq. (|5^). As before, we expect fluctuations 
to lower the scaling dimension of this coupling from its mean field value (although 
actually in d = 2 this suppression will only be logarithmic). On the other hand, 
the efficacy of the A ^ reaction is certainly not reduced by fluctuations. Hence, 



1786 P. Frdjdh, M. Howard & K. B. Lauritsen 



for rf < 2, we expect that ~ /is ~ fi^ will always run to the fixed point at oo, 
corresponding to the ordinary transition. In particular we expect this picture to 
also hold in d = 1, although in that case the transition at am = Cm, critical will not 
be accessible to controlled perturbative expansions. In d = 1 we therefore expect to 
find the O* ((T,„ < cTm, critical) Or O {am — fm, critical) transitions mentioned earlier. 
Hence we find a state of affairs very different from mean field theory: fluctuation 
effects now ensure that only the ordinary transition is accessible for d < 2, when 
f-ls > 0. 

9. Exact Results 

In the last section we reviewed the use of field theoretic methods to understand 
the effects of fiuctuations in the boundary DP and BARW models. Unfortunately, 
for the case of BARW, the presence of a second critical dimension d'^ prevents a 
controlled, perturbative investigation of the system in d = 1. Hence, given the 
fundamental difficulties associated with the BARW field theory, it seems fruitful to 
search for alternative approaches. One such alternative is provided by the theory 
of quantum spin Hamiltonians to which we now turn. 

The methods we review, first presented in Ref. are a straightforward exten- 
sion of the work in Refs. ^ and |5^. The starting point is the following set of rules 
for BARW with to = 2 in d = 1: 

074^^0 with rate 15/2 

AA 00 with rate A (56) 

0^0 ^ AAA and (!)AA ^ AA^ with rate a/2. 

Note that these rules are fermionic in character (no more than one particle per 
site is permitted) in contrast to the bosonic rules employed in the derivation of the 
earlier field theories.E3 The model described in Eq. ( p6| ) can be transformed into 
a spin picture by writing the configuration of a semi-infinite system as a vector 
\si, S2, S3, . . .), where Si = 1/2 if the i-th site is empty, and Si — —1/2 if that site is 
occupied. Hence the system ket is given by 

and the equation governing the time evolution is 

dt\p{t)) = -n\p{t)), (58) 

where, using a representation in terms of Pauli matrices, and defining = (1 — 
(T^)/2, Vk = l- Uk, = ((T^ ± ial)/2, we haveEl 

fe=i 
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+ 2X[nknk+i - 44+1]) + 2 Z^(^ " crt-i<^k+i)'^k 

fc=2 

2?i/SEP + Aff^sA + aTjBA^w (59) 



fe=l fc=l fe=2 



BARW 
fc 



Here we have used some of the notation of Ref. where SEP (symmetric exclusion 
process) refers to the diffusion piece, RSA (random-sequential adsorption) to the 
annihilation piece and BARW to the branching piece of the "quantum Hamiltonian" . 
Notice that the boundary has been included in Eq. (js^), since particles may not 
hop to the left of site 1, and the annihilation/branching processes have also been 
restricted to sites 1,2,3.... Hence, the above operator Ti, governs the evolution 
of a c? = 1 BARW system without an A ^ reaction at the boundary. Averages 
are calculated using the projection state (| = X){si}({'^*}l' ^-^^ i-^) ~ 
Following Ref. |3^, we now define an operator V where 

V = 7_i7o7i72 . • . , (60) 

with 

I2k-i = i[(l+z)a^-(l-i)], 

l2k = l[a + ^WM+l-a-^)]■ (61) 
Defining a new "quantum Hamiltonian" a.s Ti. — [D^^HV]'^ , we find 



oo 

RSA 
k 



n=iD-x)Y: /.r^^ + (o + A) E +xY^h 

k=l k=l k=l 

+ ^{ninQ - + nivo - s+Sg), (62) 

where we have used the commutation rules described in detail in Ref. Hence, 
when D = A + a, we have the following processes occurring: 

id.Ai+i ^ ASi+i rate(A + a)/2, i = 1,2,3,..., 

AiAi+i^^Si+i rate A, « = 1,2,3,..., 

(l)i^iASi+i ^ A,^iA,A,+i rate a/2, i = 1,2,3,..., 

0,_iA,A,+i ^ A,_iA,0,+i rate a/2, i = 1,2,3,..., 

AqAi 0001 rate A/2, 

00^1 Ao0i rate A/2. (63) 

Excepting the boundary terms, we see that the Hamiltonian has been mapped back 
onto itself. Furthermore, at the edge, the particles may only hop from site 1 to site 
but never the other way round. This means that we can forget about the 0-th site 
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in exchange for allowing the processes A1A2 ^ ^i02 (with rate a/2), and Ai — > 0i 
(with rate A/2). Hence, we see that the new Hamiltonian Ti corresponds to the 
case where /is 0, with the DP processes A ^ A + A and A — > generated on the 
boundary. 

If we choose the initial condition to be an uncorrelated state with density 1/2, 
denoted by |l/2}, then the density at site k, pk{t), is given by 

Pfe(i) = (|nfcexp(-m)|l/2}. (64) 

Following exactly the procedure in Rcfs. ^ and |5^ (starting with insertions of the 
identity operator VV'^ into the RHS of Eq. (H)), one can straightforwardly show 
that 

Pk{t) = i[l - (0| eM~nt)\k - 1, fc)], (65) 

where (0| is the vacuum state (with no particles), and \k ~l,k) is the initial state 
with only two particles situated at sites k — 1 and k. In Eq. ( |65| ) the LHS is the 
density at the fc-th site, whereas the RHS is one-half times the probability that a 
cluster initiated at t = by two particles at sites fc — 1 and k has not yet died out 
by time t. According to our earlier analysis, for A < 0, the LHS should scale as 
A I ''dens (far from the waU) or |A|'^i^dc„s (dose to the wall), and the RHS as \Af 
(far from the wall) or |A|'^i.='=<=<i (close to the wall). Thus, at the line D = A + a, we 
have proven the result 

(and, of course, the bulk result /Sgoed = /3dens)- We note that the bulk result was 
first proven in Ref. and a very similar result for the A + A ^ ?l reaction was 
derived in Ref. |5^ (connecting the O* and Sp* regimes). Using universality, it was 
postulated in Ref. ^ that this equality between the two surface exponents is valid 
everywhere close to the transition line, and not just where D = X + a. 

At the line D = X + a, it is now straightforward to derive a second relation 

/3l,sccd=/3Me„s- (67) 

One simply starts off with the quantum Hamiltonian TC and then follows the same 
steps as above. Ti. can then be mapped back onto the starting Hamiltonian 7i, 
meaning that the transformation is actually a duality transformation. A relation 
like that in Eq. (|6^) can then be derived, giving the above result. 

In summary, in = 1 and at the particular line in parameter space D — X + a, 
we have mapped BARW at the special transition onto BARW at the ordinary 
transition (and vice- versa), a rather non-trivial procedure reminiscent of a related 
transformation in the Ising model. Unfortunately, the results (|6^ ) and (^7|) are only 
derived for one line in parameter space and we have to rely on universality in order 
to claim that they are valid elsewhere close to the transition line. 
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10. Numerical Methods 

We now discuss the various numerical methods which can be used to obtain precise 
estimates for the critical exponents for boundary DP and BARW. The values of the 



exponents are listed in Section 10.3 



10.1. Extraction of exponents from Monte- Carlo simulations 

In this subsection we review how the exponents are extracted from Monte-GaiJn^im= 
ulation data for the DP and DP2 models with walls in d = 1 and d = 2Jil't3fiilO 
Numerics have also been carried out for BARW with Eiwall in d = 1, with re- 
sults in good agreement with boundary DP2 simulationsEJ The most efficient way 
of extracting information on boundary DP and DP2 is to perform simulations at 
criticality, starting from an inital seed next to the wall. Measurements can then be 
taken of the survival probability Pi (t) , the activity in the bulk Ni (t) and at the wall 
Ni_i{t), the average spread, of the cluster and the probability pi{s) to have 

a cluster of size (mass) sSB Furthermore, by averaging over surviving clusters only 
(denoted by an over-line), measurements can be taken of the surviving bulk activity 
Ni{t) and the surviving wall activity iVi.i(t). Although carrying out these simula- 
tions is straightforward, extracting the exponents discussed in the previous sections 
requires some further analysis. _We will now review the scaling theory which allows 
these connections to be madeO Note that all the relations given below are valid 
for the DP ordinary transition and for both the IBC (special) and RBC (ordinary) 
DP2 transitions, and hence these labels will be suppressed from now on. 

The probability for a cluster grown from a seed on the wall still to be alive at 
time t is given by Eqs. ( p^ and (^^. At criticality (A = 0) it has the following 
behavior 

Pl{t) = t-*l--V(tAc), <5l,SCCd = /3l,SCCd/i^||, (68) 

where the scaling function / depends on the cutoff tc ^ C|| ^ lAj""^'!. Hence, 
the probability of growing a cluster which lives exactly t time steps behaves as 
pi{t) ~ t-i-'Si.-cd^ 

The average number of active sites at criticality, averaged over all clusters, is 
obtained by integrating the density ( p7| ) and ( |38| ) over space, giving 

Niit)^^', Kl dx - ^dcns - <^l.sood, (69) 

where we have introduced the cluster envelope or "roughness" exponent x — ^'±/i'\\ 
(= 1/z), and the notation Jdcns — Pdcns/i^\\- Note that the scaling relation in 
Eq. ( |69| ) corresponds to the hyperscaling relations (^0|) and (^), a fact which follows 
from {si{t)) = dt' Ni{t'), and the relation 71 — t^||(l -t- ki). By integrating the 
density on the wall (^8|) and (^9|), we obtain the average number of active sites at 
criticality on the wall 

Ni^lit) ^ t^'-\ Ki,i = {d-l)x~ (5l,do„s - ^l,seed, (70) 
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where ^i,dons — /3i,dcns/i^|| • Note also that the scahng relation in Eq. ( [70| ) corre- 
sponds to the hyperscaling relations ( ^3| ) and (^3|) at criticality, since (si,i(t)) = 
/odt'iVi,i(t'), and 7i_i =z/||(l + «:i,i). 

Alternatively, by averaging only over clusters which survive to infinity (denoted 
by an over-line), we obtain 

'Thit) ^t~, Kl^dx- <5dc„s. (71) 

The activity on the wall for surviving clusters reads 



iVi,i(i) kTJ = (d - l)x - <5i,dcns. (72) 

Simulations in d = 1 thus directly yield ^i^dcns- The average position of activity 
follows from Eqs. ( p7| ) and (38), giving (x^) ~ i^^, where x is the distance from the 
seed and the average is taken over all active points at a given time. 

Additional numerical data was obtained in Refs. 21 2^ by considering the cluster 
size distributions at criticality. In the bulk the typical cluster size Sc of finite clusters 
scales as volume times density, i.e., 

Sc^CU\MA)^\A\-'/^, l/a = diy^+i^\\-pdens- (73) 

From the lifetime survival distribution (|68|), it is then straightforward to obtain 
the probability Pi{s) to have a cluster of size larger than s, for clusters started 
from a seed on the wall. Using the fact that the lifetime is set by the parallel 
correlation length, t ^ lAj""^!! , we see that the typical cluster size and lifetime 
are connected by s ~ i^/'^ii'^. Hence we obtain Pi{s) ^ Pi{t - s''"'^) ~ s-/3i,=codO-^ 
Thus, we eventually obtain the probability pi(s) to have a cluster of exactly size s, 
Pi{s) — —dPi{s)/ds, with the result 

Pi{s) = s-^^g{s/s,), m = l+ , ■ (74) 

dl^± + I^W - /5dcns 

Similarly, the cluster size distribution on the wall due to a seed located at the wall 
can also be obtained. In this case the typical cluster size of finite clusters is 

Swaii,c~^f'^||ni(A)^|Ari/'^S (75) 

where the cutoff exponent is 

1/fTi = (rf - l)i^_L + Z^ll - /?l,dcns- (76) 

The resulting distribution reads 

Pl,l(Swau) = S^rn '/(Swall/Swall.c), (77) 

with 

1 I /?l,SCCd /'70\ 

[d - + i^ll - /^i,dcns 
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With the above scahng relations it is now straightforward to extract all the 
exponents from the numerical data. Values for the exponents obtained in this way 
are listed in Section 10. 3| . An additional interesting exponent relation for the RBC 



transition for DP2 can be obtained by assuming that the survival probability is 
dominated by the return to the wall of the cluster-envelope. This leads toliS 



«d = l-X, (79) 



in agreement with the simulation results for the RBC (see Section 10. 3| ). Quail 



tatively, this means that the I2 regions located at the wall determine the scaling 
since they can only disappear when the activity returns to the wall. Note that a 
relation of this kind clearly fails for the IBC transition. Furthermore, if the cluster 
lifetime is defined to be the return time of the cluster-envelope (i.e. the return time 
of the rightmost active site) to the initial point, then we expect clusters defined in 
this way to have a lifetime distribution exponent (5i,scod given by Eq. ([79|). This 
prediction is in agreement with the simulations in Ref. where various models in 
the DP and BARW classes were studied with cluster lifetimes defined in the way de- 
scribed above. For further information on the DP2 boundary exponents, including 
inequalities for the /3i,sGed and /3i,dons exponents, we refer to Ref. 

10.2. Series expansions 

In the previous subsection we reviewed the scaling relations used to extract ex- 
ponents from Monte-Carlo simulations. These simulations provide fairly accurate 
exponent estimates for both the DP and BARW classes (see the next subsection). 



However, for the case of DP in d = 1 (but curiously not for BARW in d = Icj) this 
level of accuracy can be far surpassed by using the method of series expansions JH §3 
With uncertainties in the seventh or eighth digit, these exponent values are very 
reliable and useful for testing conjectures and scaling laws. For example, most of 
the DP exponents in d = 1 are now known not to lie close to any simple rational 
numbers. In fact, it seems more likely that the exponents are irrational. However, 
one important exception is the exponent n governing the mean cluster lifetime in 
the presence of a wall (see Eq. (|25|)). This exponent has been conjectured to equal 
unityjllalthough this has now been challenged by the estimate ri = 1.00014(2) (see 
below) E3 

The idea behind series expansions in DP is to find an effective algorithm for 
calculating the expansion coefficients of moments of the pair-connectedness (the 
probability that the site x at time t is occupied given a seed at the origin at i = 0). 
By analyzing these expansions, critical parameters of the model can be obtained, 
including estimates for pc, and the exponents 7, i^^, and . These series expansions 
can be generated by transfer matrices that relate states at time t to those at one 
time step later. At any time t this gives polynomials in the percolation probability p. 
Since the difficulty in generating these series expansions is of exponential complexity, 
a large amount of effort has been devoted to keeping the computational effort as 
small as possible (see Ref. ^ and references therein). Once a series expansion is 
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derived, it is then analyzed by using differential approximants. The more terms that 
can be computed in the series, the more accurately the exponents can be estimated. 

10.3. Exponent Values 

In this section we list the best exponent values currently available for the surface 
and bulk exponents for DP and BARW. As we have mentioned, for DP in d = 1 
the best results are from series expansions; in all other cases we give Monte-Carlo 
data. Note that various attempts have been made to fit these exponent values with 
rational numbers; further details can be found in Refs. p^ , |27| . 

Table 4. Critical exponents for DP in d = 1 and d = 2 in the bulk and for the ordinary transition 
at the boundary. 





d = l 


d = 2 


Mean Field 


f^dcns — 




0.159 464(6) 


0.450(1) 


1 


f^l ,dens 


*^l,secd 


0.423 17(2) 


0.82(4) 


3/2 


K 




0.313 686(8) 


0.230(1) 





Kl 




0.049 98(2) 


-0.13(1) 


-1/2 






0.473 150(10) 


0.682(1) 


1 






1.108 247(10) 


1.268(8) 


3/2 


Ml 




1.287 25(2) 


1.49(8) 


7/4 


2X 




1.265 226(13) 


1.132(7) 


1 


0dens — 




0.276 486(8) 


0.583(4) 


1 


01, dens 


/5l,sccd 


0.733 71(2) 


1.07(5) 


3/2 


T 




1.457 362(14) 


0.711(7) 





n 




1.000 14(2) 


0.26(2) 





7 




2.277 730(5) 


1.593(7) 


1 


71 




1.820 51(1) 


1.05(2) 


1/2 



Table 5. Critical exponents for the BARW class measured from simulations of DP2 in d = 1. 





d = l 


d=l (IBC) 


d = l (RBC) 




0.287(5) 


0.288(2) 


0.291(4) 


*^sccd 


0.290(5) 






'5l,sGcd 




0.641(2) 


0.426(3) 


'^1 , dens 




0.415(3) 


0.635(2) 


K 


0.000(2) 






Kl 




-0.354(2) 


-0.141(2) 


= "ki 


0.288(5) 


0.287(2) 


0.285(2) 




1.225(5) 






Ml 




1.500(3) 


1.336(3) 


Ml,i 


1.408(5) 


2.05(5) 


2.15(5) 


2X 


1.150(5) 


1.150(3) 


1.152(3) 


Pdcns 


0.922(5) 


0.93(1) 


0.94(2) 


/-^sccd 


0.93(5) 






/-^l .seed 




2.06(2) 


1.37(2) 


/^l ,dens 




1.34(2) 


2.04(2) 


T 


2.30(3) 






ri 




1.16(4) 


1.85(4) 


7 


3.22(5) 






71 




2.08(4) 


2.77(4) 


71,1 


1.38(3) 


« 0) 


« 0) 
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11. Other Directions 

In this section we briefly mention some other nonequilibrium systems where bound- 
ary effects have been analyzed. Additional studies, which we w ill n ot discuss, include 
(i) the connection between surface DP and local persistencepj (ii) d = 1 density 
matrix renormalization group calculations of some reaction-diffusion processes;B3 

(iii) active, but slanted, walls in DP, which give rise to a "curtain" of activity 
whose width is given by an angle-dependent correction to bulk DPpj and, finally, 

(iv) Monte-Carlo simulation studies of rigidity percolation with and without walls, 
which have shown the model to belong to the DP universality class.El 



11.1. Compact DP 

By forcing a DP cluster to be dense one obtains a different universality class called 
Compact Directed PercolationEj This actually corresponds to the special case q = 1 
of the Domany-Kinzel model in Figures 0and I The constraint q = I forces all 
interior sites and/or bonds to be active, whereas the size of the cluster is controlled 
by p, which now only affects the boundaries of the cluster. As the dynamics of the 
cluster is entirely controlled by the dynamics of its surface, the problem is greatly 
simplified. For instance, in c? = 1, the extinction of a compact DP cluster can be 
viewed as the probability for a pair of random walkers to coincide. It is therefore no 
surprise that compact DP can be exactly solved in c? = 1 and that all the exponents 
are simple integers.l23 For bulk compact DP in d = 1, one has /3scod = 1, = 2, 
h'^ = 1, and 7 = 2. 

By introducing a wall in compact DP, the survival probability is altered and one 
obtains surface critical exponents just as for DP. With an inactive wall, the cluster 
is free to approach and leave the wall, but not crosa^Fpr d = I, this gives rise 
to /3i,seed = 2, i.e., twice as big as /3sccd for the bulk.E3'Ea On the other hand, for 
an active wall, the cluster is stuck to the wall and therefore described by a single 
random walker for d — I. By reflection in the wall, this may be viewed as symmetnc 
compact DP which has the same /3soGd as normal compact DP, giving /Ji^seod = lEI 



11.2. Backbone and red bonds in DP 

In Ref. [l^ the so-called backbone and red bonds of d = 1 DP clusters have been 
investigated. The backbone is obtained from the infinite cluster by removing all 
dangling ends. Thus the backbone consists of precisely those bonds which would be 
occupied by both the time-directed DP process and its reverse time-directed process. 
It then follows that the backbone density |A|'^ is described by the exponent 
[3^^ = 2/3. Numerically, the backbone dimension on the wall was measured with 
the result Df^ = 0.16 ±0.01. Using the scaling relation Df ^ = 1 -/3f it then 
follows that /3f^ = 1.46 in good agreement with /3f^ — 2/3i, cf. the result for the 
bulk. 

On the backbone one can identify so-called red bonds0 that^jf one is cut, divide 
the cluster into two parts. A renormalization group argumentEj (see also Ref. |63|) 
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yields that the number of red bonds up to time t scales as 

NR{t) ^ t^/^^ii . (80) 

Reference ^ measured the scaling of red bonds for DP with a wall and obtained 
results in complete agreement with Eq. (|^. In addition the scaling of red bonds 
along a longitudinal cut for DP with no wall was measured with the result ./VJj"* (t) ~ 
^-o.o4±o.02^ This is in accordance with the expected behavior 7V];,"*(t) ^ t^/''\i-''^/''ii ~ 
^-0.056^ where the extra factor originates from the scaling of the width ^ t'^^/'^ii 
of the cluster. Finally, the scaling of boundary red bonds for DP with the wall was 
measured, where the behavior Nii^i{t) ~ ^-o.eoio.i ^g^g found. 

11.3. Invasion Percolation 

Finally, we briefly discuss a slightly different system, namely invasion percolation 
(IP), which is a model for the growth of infinite percolation clusters at criticalityO 
In d = 2 the fractal dimension of IP clusters is D = 2 - = 2 - (5/36)/(4/3) = 
91/48 ~ 1.90, where the known exact values for percolation in d = 2 have been used 
(see, e.g., Ref. p5| ). 

In Ref. ^6|, boundary effects in the growth of d = 2 invasion percolation clusters 
were studied. Numerically, near a wall it was found that the fractal dimension of IP 
clusters was Ds = 1.67±0.03. Using the scaling theory for boundary nonequilibrium 
systems, it follows that the activity on the wall has the fractal dimension Di = 
1 — (3i/i^, where /?i is the surface exponent for d = 2 percolation. Using f3i — 4/9, 
we obtain that Di — 2/3. The fractal dimension of parts of the cluster close to the 
wall is thus predicted to be Di + 1 = 5/3 in excellent agreement with the simulation 
result. 

12. Summary 

In this review we have outlined the boundary critical behavior of some nonequilib- 
rium systems, with a particular focus on the directed percolation and branching- 
annihilating random walk universality classes. Through the use of a variety of 
theoretical and numerical techniques, including mean field, scaling and field theo- 
ries, exact solutions, Monte-Carlo simulations, and series expansions, a considerable 
amount of progress has been made in understanding the boundary critical proper- 
ties of these models. Nevertheless some important problems remain open and we 
wish to conclude this review with a brief list of some of these remaining questions. 

1. Due to the presence of a second critical dimension d[, the BARW field theory in 
d = 1 is not well-controlled either at the boundary or in the bulk. An improved the- 
ory would be highly desirable in confirming the results provided by exact solutions 
and Monte-Carlo simulations. 

2. Although we have not discussed it in this review, the boundary critical prop- 
erties of dynamical percolation form another interesting case (for d > 2).E3 In this 
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universality class the development of an appropriate boundary field theory, and the 
calculation of boundary exponents using epsilon expansion techniques, has not yet 
been attempted. 

3. The intriguing question of whether ti = 1 for the DP ordinary transition in = 1 
remains unsolved. No theoretical explanation has emerged for why ti should equal 
unity, and, in fact, the latest series expansions yield a value for n very slightly away 
from Ti = 1. Hence the possibility of numerical coincidence, where ri just happens 
to lie very close to unity, cannot be ruled out. 

4. Most of the numerical work so far has focused on the ordinary transition. Par- 
ticularly for the case of DP, it would be useful to measure the exponents at the 
extraordinary and special transitions for d > 1. 

Acknowledgements 

Part of this work was carried out while M.H. and K.B.L. were at the Niels Bohr 
Institute. P.F. acknowledges support from the Swedish Natural Science Research 
Council. M.H. acknowledges support from the CATS group at the Niels Bohr Insti- 
tute, from the NSF through the Division of Materials Research and from NSERC 
of Canada. K.B.L. acknowledges support from the Carlsberg Foundation. 

References 

1. K. Binder, in Phase Transitions and Critical Phenomena, Vol. 8, edited by C. Domb 
and J. L. Lebowitz (Academic Press, London, 1983). 

2. H. W. Dichl, in Phase Transitions and Critical Phenomena, Vol. 10, edited by C. 
Domb and J. L. Lebowitz (Academic Press, London, 1986). 

3. H. W. Diehl, Int. J. Mod. Phys. B 11, 3593 (1997). 

4. W. Kinzel, in Percolation Structures and Processes, edited by G. Deutscher, R. 
Zallen, and J. Adler, Annals of the Israel Physical Society, Vol. 5 (Adam Hilger, Bristol, 
1983). 

5. R. Dickman, in Nonequilibrium Statistical Mechanics in One Dimension, ed. V. 
Privman (Cambridge University Press, Cambridge, 1997). 

6. H. Hinrichsen, Adv. Phys. 49, 815 (2000). 

7. H. Takayasu and A. Yu. Tretyakov, Phys. Rev. Lett. 68, 3060 (1992). 

8. I. Jensen, Phys. Rev. E 50, 3623 (1994). 

9. J. Cardy and U. C. Taubcr, Phys. Rev. Lett. 77, 4780 (1996); J. Stat. Phys. 90, 1 
(1998). 

10. P. Grassberger, F. Krause, and T. von der Twer, J. Phys. A: Math. Gen. 17, L105 
(1984); P. Grassberger, J. Phys. A: Math. Gen. 22, L1103 (1989). 

11. M. H. Kim and H. Park, Phys. Rev. Lett. 73, 2579 (1994); H. Park, M. H. Kim, and 
H. Park, Phys. Rev. E 52, 5664 (1995). 

12. W. Hwang, S. Kwon, H. Park, and H. Park, Phys. Rev. E 57, 6438 (1998). 

13. W. Hwang and H. Park, Phys. Rev. E 59, 4683 (1999). 

14. N. Menyhard and G. Odor, J. Phys. A 29, 7739 (1996). 

15. H. Hinrichsen, Phys. Rev. E 55, 219 (1997). 

16. H. K. Janssen, B. Schaub, and B. Schmittmann, Z. Phys. B 72, 111 (1988). 



1796 P. Frdjdh, M. Howard & K. B. Lauritsen 



17. J. W. Essam, A. J. Guttmann, I. Jensen, and D. TanlaKishani, J. Phys. A 29, 1619 

(1996) . 

18. K. B. Lauritsen, K. Sneppen, M. Markosova, and M. H. Jensen, Physica A 247, 1 

(1997) . 

19. P. Frojdh, M. Howard, and K. B. Lauritsen, J. Phys. A 31, 2311 (1998). 

20. H. Hinrichsen and H. M. Koduvely, Eur. Phys. J. B 5, 257 (1998). 

21. K. B. Lauritsen, P. Frojdh, and M. Howard, Phys. Rev. Lett. 81, 2104 (1998). 

22. L Jensen, J. Phys. A 32, 6055 (1999). 

23. M. A. de Menezes and C. F. Moukarzel, Phys. Rev. E 60, 5699 (1999). 

24. E. Carlon, M. Henkel, and U. SchoUwock, Eur. Phys. J. B 12, 99 (1999). 

25. C.-C. Chen, H. Park, and M. den Nijs, Phys. Rev. E 60, 2496 (1999). 

26. M. J. E. Richardson and Y. Kafri, Phys. Rev. E 59, R4725 (1999); Y. Kafri and M. J. 
E. Richardson, J. Phys. A 32, 3253 (1999). 

27. M. Howard, P. Frojdh, and K. B. Lauritsen, Phys. Rev. E 61, 167 (2000). 

28. P. Grassberger, J. Phys. A 25, 5867 (1992). 

29. J. L. Cardy and R. L. Sugar, J. Phys. A 13, L423 (1980). 

30. H. Janssen, Z. Phys. B 42, 151 (1981). 

31. P. Grassberger and K. Sundermeyer, Phys. Lett. B 77, 220 (1978). 

32. E. Domany and W. Kinzel, Phys. Rev. Lett. 53, 311 (1984). 

33. W. Kinzel, Z. Phys. B 58, 229 (185). 

34. Note that the convention used for A here is different from that in Ref. ^ 

35. P. Grassberger and A. de la Torre, Ann. Phys. (N.Y.) 122, 373 (1979). 

36. K. Mussawisade, J. E. Santos, and G. M. Schiitz, J. Phys. A 31, 4381 (1998). 

37. L Jensen, J. Phys. A 29, 7013 (1996). 

38. J. W. Essam, K. De'BeU, J. Adler, and F. M. Bhatti, Phys. Rev. B 33, 1982 (1986). 

39. J. W. Essam, A. J. Guttmann, and K. De'BeU, J. Phys. A 21, 3815 (1988). 

40. P. Grassberger, J. Phys. A 22, 3673 (1989). 

41. P. Grassberger and Y.-C. Zhang, Physica A 224, 169 (1996). 

42. L Jensen, Phys. Rev. A 45, R563 (1992). 

43. M. A. Mufioz, R. Dickman, A. Vespignani, and S. Zapperi, Phys. Rev. E 59, 6175 
(1999). 

44. L Jensen, J. Phys. A 32, 5233 (1999). 

45. L Jensen, J. Phys. A 30, 8471 (1997). 

46. D. Zhong and D. ben-Avraham, Phys. Lett. A 209, 333 (1995). 

47. B. M. McCoy and T. T. Wu, Phys. Rev. 162, 436 (1967). 

48. A. J. Bray and M. A. Moore, J. Phys. A. 10, 1927 (1977); T. W. Burkhardt and H. 
W. Diehl, Phys. Rev. B 50, 3894 (1994). 

49. J. L. Cardy, Nucl. Phys. B 240 [FS12], 514 (1984). 

50. H. A. Kramers and G. H. Wannier, Phys. Rev. 60, 252 (1941); ibid. 60, 263 (1941). 

51. A. Drewitz, R. LeidI, T. W. Burkhardt, and H. W. Diehl, Phys. Rev. Lett. 78, 1090 
(1997). 

52. Note that the mean field surface density exponent /3j^^,js quoted in footnote [22] of 
Ref. ^ is incorrect. 

53. J. L. Cardy, J. Phys. A 16, 3617 (1983). 

54. M. A. Munoz, G. Grinstein, and Y. Tu, Phys. Rev. E 56, 5101 (1997). 

55. L. Peliti, J. Physique 46, 1469 (1985). 

56. G. M. Schiitz, Z. Phys. B 104, 583 (1997). 

57. Note that some of our exponent definitions in this section differ from those sometimes 
found in the literature. In particular our k exponents are sometimes called 77, and the 
dynamic exponent z is sometimes redefined as z — > 2/z. 

58. R. Bidaux and V. Privman, J. Phys. A 24, L839 (1991). 



Directed Percolation and Systems with Absorbing States 1797 



59. J.-C. Lin, Phys. Rev. A 45, R3394 (1992). 

60. J. W. Essam and D. TanlaKishani, J. Phys. A 27, 3743 (1994). 

61. H. E. Stanley, J. Phys. A 10, L211 (1977); R. Pike and H. E. Stanley, ibid. 14, L169 
(1981). 

62. A. Cogniho, Phys. Rev. Lett. 46, 250 (1981); J. Phys A 15, 3824 (1982). 

63. G. Hubcr, M. H. .Jensen, and K. Sneppen, Phys. Rev. E 52, R2133 (1995). 

64. T. Vicsek, Fractal Growth Phenomena, 2nd ed. (World Scientific, Singapore, 1992). 

65. M. Henkel, Conformal Invariance and Critical Phenomena (Springer- Verlag, BerUn, 
1999). 

66. R. Cafiero, G. CaldarelU, and A. GabrieUi, Phys. Rev. E 56, R1291 (1997). 



